Long-range repulsive interaction between TTF molecules on a metal surface 

induced by charge transfer 
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The low-coverage adsorption of a molecular electron donor, tetrathiafulvalene, on Au(lll) is 
characterized by the spontaneous formation of superlattice of monomers, whose spacing exceeds 
the equilibrium distance of non-covalent interactions and depends on coverage. The origin of this 
peculiar growth mode is due to a long-range repulsive interaction between molecules. The analysis 
of molecular-pair distributions obtained by scanning tunneling microscopy measurements permits 
us to determine that the nature of TTF intermolecular interactions on Au (111) is electrostatic. A 
repulsion between molecules is caused by the accumulation of charge due to electron donation into 
the metal surface, as pictured through density functional theory calculations. 



The spontaneous formation of self-organized molecular 
structures at metal surfaces follows a complex balance of 
interactions between the basic functional units [1 . At- 
tractive short-range forces between molecules are ubiq- 
uitous during growth, but their strength and relevance 
varies depending on the molecular functionalization. For 
the case of adsorption on metal surfaces these forces com- 
pete with substrate- mediated interactions, for example, 
through elastic stress fields [2j [3l [4] or through surface 
state electrons O IH [3 [H [H [10] . These usually extend for 
larger length scales than intermolecular dispersion forces 
and can lead to characteristic quasi-periodic arrays of 
particles [H [9] . Long-range interactions can also have a 
repulsive nature. This is the case of electrostatic interac- 
tions between charged particles weakly interacting with a 
non-conducting host support [HI [13] , or in ensembles of 
organic molecules with large dipolar moments on metal 
surfaces [T5]. 

Apolar and neutral molecules are not expected to build 
up long-range interaction potentials other than those me- 
diated by the underlying substrate and, in most 
cases, attractive dispersion forces lead to nucleation in 
two or three dimensional condensates. Charge redistri- 
bution upon molecular chemisorption is also able of ren- 
dering interesting changes in the interaction potentials 
between molecules [16]. Although this effect is presum- 
ably strong in charge transfer adsorbate systems, it has 
been usually neglected due to the screening nature of 
metallic substrates. An experimental proof of its rele- 
vance in intermolecular interactions is thus still missing. 
This could also help to build up a quantitative picture 
about fundamental processes related to molecular charg- 
ing on metal surfaces. 

Here, we report the spontaneous formation of 
quasiperiodic superlattices of single tetrathiafulvalene 
(TTF) molecules on a Au(lll) surface driven by lo- 
cal charges at the interface induced upon chemisorption. 



TTF is well known as a prototype donor molecule in 
charge transfer compounds [17 . The free molecule has 
no electrical dipole moment. However, on Au(lll), it be- 
comes charged upon electron donation. Using a combina- 
tion of low-temperature scanning tunneling microscopy 
and density functional theory (DFT) we resolve that a re- 
pulsive long-range interaction between charged molecules 
is built up, thus hindering nucleation in islands. Through 
the analysis of molecular pair distributions we recon- 
struct the coverage-dependent intermolecular potential 
wells forming the molecular lattice. 

The experiments were carried out in a custom-made 
low temperature STM under ultra-high vacuum. An 
atomically clean Au(lll) substrate is exposed to a con- 
tinuous flux of TTF molecules sublimated from a home- 
made Knudsen cell. TTF has a very high vapor pressure. 
To obtain the low coverages used here the crucible was 
heated to 30° C. The sample temperature was varied dur- 
ing the dosing between 80 K and 300 K, and posteriorly 
cooled down to the operating temperature of the STM 
(4.8 K) for sample inspection. 

Room temperature deposition of a small amount of 
TTF (< 0.1 ML) leads to the formation of a characteristic 
quasi-periodic one-dimensional array of TTF monomers 
along the FCC regions of the Au(lll) 23x^/3 recon- 
struction (Fig. [1]). The separation between monomers 
amounts several nanometers, ~ 3 nm for the data in 
Fig. [ija). This distance is significantly larger than 
the typical length scale of attractive non-covalent in- 
teractions. The formation of the superlattice of TTF 
monomers needs to be thermally activated. Fig. 1(b) 
shows the result of dosing TTF on a 80 K cold sample. 
In this case both monomers and small TTF clusters ap- 
pear randomly spread and are easily dragged by the STM 
tip, probably because they populate a weakly adsorbed 
precursor state. Only upon annealing the molecules self- 
organize forming the distinctive quasi-periodic array, as 



FIG. 1: (a) STM image of an Au(lll) region with 0.03 ML 
of TTF deposited at room temperature. One-dimensional ar- 
rays of TTF monomers follow the Au(lll) herringbone re- 
construction . (b) Depositing molecules on a cold sample (80 
K) leads to population of a weakly adsorbed precursor state, 
in which molecules may nucleate in clusters. After annealing 
to room temperature (c), the TTF arrays along FCC regions 
are formed, (d) STM image (inset; Vs=-1 V) and its Laplace 
filtered image 30 of a TTF molecule. The latter reveals that 
two of the sulfur atoms appear brighter suggesting a small 
tilt of the molecular plane with respect to the surface (later 
confirmed by theory). 

shown in Fig. [TJc). In this case, high-resolution STM 
images (Fig.jljd)) of intramolecular structure can be ob- 
tained. At negative sample bias we find TTF monomers 
as composed by four protrusions. The two largest corre- 
spond to the 4 sulfur atoms. The other two are fainter 
tails due to the ethylene ends. TTF appears with a char- 
acteristic asymmetry in the images, resembling two of the 
lateral S atoms being higher than the other two. 

Our results clearly indicate that, in the low-coverage 
limit and after annealing, TTF does not respond to at- 
tractive forces like, for example, hydrogen bonding to 
sulfur atoms [18], avoiding nucleation into islands. Such 
behavior prevails as the coverage is increased, accom- 
panied by a monotonous decrease in the average pair 
distance (Fig. [2|a-c)). At 0.08 ML the array is com- 
pressed (average pair distance ~ 2nm) into double rows 
of monomers in the FCC regions of the reconstruction. 
Close to this coverage HCP regions start also to be 
populated with similar one-dimensional arrays of TTF 
monomers. Such tendency to avoid nucleation through 
the formation of quasi-periodic molecular arrays is in- 
dicative of a long-range interaction mechanism different 
from (shorter-range) non-covalent dispersion forces be- 
tween molecules. 

Elastic deformation of the substrate can lead to long- 
range interactions between adsorbates [21 [3l [5] . The in- 



FIG. 2: (a-c) STM images of TTF on Au(lll) at various 
coverages. At 0.08 ML molecules appear already at the HCP 
regions, (d) Pair distributions / of the one-dimensional TTF 
arrays for the data shown in Fig. 1 (0.03 ML) and Fig. 2(a- 
c). For 0.08 ML and 0.16 ML the distributions are performed 
on HCP regions. More than 500 pairs are analyzed in each 
plot. The molecular coverage is determined from STM images 
of large surface areas, assuming that 1 ML corresponds to 2 
molecules/nm^. From the lowest to the largest coverage we 
obtain an average pair distance r of 3.5 nm, 2.5 nm, 3.3 nm, 
and 1.7 nm in the one dimensional arrays. The corresponding 
ID distribution functions for non interacting particles fran 
are included. 



duced stress field can oppose the approach of two ad- 
sorbates becoming the driving force of an ordered phase. 
Indeed, the Au (111) herringbone reconstruction is itself 
a stressed atomic layer, and therefore, the periodicity 
of its folding is very sensitive to small changes in the 
stored elastic energy. Our data show that the herring- 
bone structure is unaffected by a sub-monolayer cover- 
age of TTF. Therefore, this mechanism is improbable in 
our case. An alternative mechanism for long-range inter- 
action between atoms [7l |8] or molecules [10] on metal 
surfaces is the interaction potential associated with the 
Friedel oscillations due to the scattering of surface-state 
electrons with the adsorbates. A key element in this 
mechanism is the oscillatory character of the interaction 
with a period related to half of the electronic Fermi wave- 
length (Af/2). For Au(lll) this corresponds to 1.8 nm, 
much smaller than the average pair distance of the data 
in Fig. [ij Furthermore, the average pair distance de- 
creases monotonously with increasing TTF density along 
the rows (Fig. [2|. Thus, an interaction mediated by sur- 
face electrons can also be discarded as the driving force 
leading to the superlattice formation. 

Fig. [2|d) shows the pair distance r distributions of one- 
dimensional arrays for various coverages (along FCC re- 
gions or HCP regions depending on the coverage). For a 
one-dimensional system of non-interacting particles the 
first-neighbors' random pair distribution function fran 
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FIG. 3: Results from DFT simulations, (a) Fully relaxed 
configuration of TTF on Au(lll). The uppermost two gold 
layers as well as the molecular degrees of freedom are relaxed 
until atomic forces are lower than 0.01 eV/A. (b) Tersoff- 
Hamman constant current image [23] of the molecule in (a) 
(V=-0.5 V). (c) Induced electronic density by the molecule- 
surface interaction, (d) Lateral (x-y planes) integration of the 
induced charge. The arrows show the vertical distance values 
at which the two topmost surface layers and the two binding 
S atoms lie. (e) Accumulated induced dipole. Together with 
(d) it reveals that the molecule becomes positively charged, 
(f) Projected density of states on molecular orbit als. The 
electronic states with HOMO character are partially empty, 
in agreement with the data of (c-e). 



decays monotonically with the pair distance r as shown 
by a dashed line in figure [19 . The peaked distributions 
in Fig. [2|d) are symptomatic of a repulsive long range 
interaction between monomers. Motivated by the donor 
nature of TTF as a free molecule we have performed ab- 
initio calculations in order to trace back the nature of 
the molecule-surface interaction and its effect in the long 
range repulsion between TTF monomers. 

We have used density functional theory (DFT) within 
the generalized gradient approximation [20^ as imple- 
mented in the VASP code [21 to evaluate the properties 
of a relaxed layer of TTF on an artificial FCC (111) 4- 
layers slab of gold atoms. The electron-ion interaction 
is described by the projector- augmented wave (PAW) 
scheme [22 . Figure |3j a) shows the resulting structure 
of a relaxed TTF molecule in a 6x4 unit cell. This large 
unit cell is employed in order to account for large molecu- 
lar separations within computationally reasonable limits. 

The interaction of the molecule with the surface is 
driven by local S-Au bonds. Due to the incommensurate 
dimensions of molecule and surface, the local interactions 
lead to an asymmetric chemisorption of the molecule. As 



a result, the molecule aligns along the {110} direction 
of the surface and tilts 8° with respect to the surface 
plane. The tilt is responsible of the asymmetry in the 
experimental constant current STM image. Fig. [TJd), as 
it is here captured by its Tersoff-Hamman simulation [23] 
(Fig. [sj^b)). At negative bias voltage, the STM image is 
basically dominated by the shape of the HOMO. 

The local interaction character between Au and S 
atoms implies a sizable bonding strength and a large 
charge donation into the surface. Indeed, the adsorp- 
tion energy after dipole corrections is -0.86 eV, and the 
surface-molecule distance is 2.76 A. The electronic struc- 
ture of the S atoms has a large contribution in the highest 
occupied molecular orbital (HOMO), what causes a large 
redistribution of electronic charge (Fig.jsjc)). The charge 
donation is expressed by a partial decrease of the elec- 
tron density in the whole molecular plane. The result 
is a positive charging of the molecule and the creation 
of a surplus of negative charge localized close to the S- 
Au bonds. Figure |3|d) shows the planar integration of 
charge. An excess of positive charge (~ 0.6 e) is located 
about the molecule and the corresponding screening neg- 
ative charge (~ —0.4 e) is between the molecule and the 
first atomic layer. The molecule-surface interaction leads 
to a large surface dipole that is evaluated in Fig. |3je) ac- 
cording to Ref . [24 . The dipole is zero inside the surface 
and builds up across the molecule reaching a value of 5.0 
D. 

The charge donation gleaned from the induced elec- 
tronic density causes the partial emptying of the HOMO. 
This is clearly seen by plotting the projection of 
the full electronic structure onto the molecular or- 
bitals corresponding to the present molecular conforma- 
tion (Fig. [sjf)) [25]. The molecule-surface interaction 
also broadens the molecular features associated to the 
HOMO-1, HOMO and the lowest unoccupied molecu- 
lar orbital (LUMO), revealing a substantial hybridiza- 
tion with the surface electronic structure, while higher- 
lying orbitals are thinner showing their small role in the 
molecular-surface interaction. 

The ab-initio results evidence a significative charging 
of the TTF on the Au(lll) surface. For pair distances 
shorter than the Thomas- Fermi screening length on gold 
surfaces this local charge can build up a repulsive poten- 
tial between molecules p6^. To evidence its role in the 
formation of the arrays we analyze the statistics shown 
in Fig. 2(d). The experimental pair distributions / arise 
from the site occupation as dictated by the Boltzmann 
factor exp(-(a;(r) — /i)//c^T), where u;{r) is the mean in- 
teraction potential behind the formation of the superlat- 
tice, /i is a (coverage dependent) zeroth order potential 
[27] and ks the Boltzmann constant. To evaluate uj{r) 
we divide the experimental pair distribution / by that 
of non-interacting particles (fran) and plot -\n{f / fran) 
(Fig. [4|. In the limit of a very dilute system [71 181 [28] . 
i.e. where no quasiperiodic array is formed, cj(r) would 
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FIG. 4: Mean interaction potentials uj{r) of one- dimensional 
TTF arrays obtained from the pair distributions shown in 
Fig. [ijd). The dashed line represents the pair electrostatic 
interaction E(r) between particles charged with 0.3 e and a 
temperature (T=160 K) to fit the repulsive part of a;(r) for 
the most dilute case. Each curve has been shifted upwards an 
amount (8.4, 5.4, 4.1, 3.8, from top to bottom) representing 
the coverage dependent zeroth order internal potential and 
approximated here as the electrostatic energy per molecule in 
a fully periodic lattice and using the fitted temperature, for 
consistency. 



be a good approximation to the (repulsive) pair interac- 
tion potential, E(r). Here, however, uj{r) has the shape of 
a potential well. As the molecular density increases the 
well becomes more symmetric and shallower, in accord 
with the TTF molecules being confined into sharper pair 
distribution and, hence, forming a superlattice. Unfor- 
tunately, it is not trivial to obtain the shape of the pair 
interaction E(r) from the mean potential Lj{r) [29 . How- 
ever, we note that for small pair distances a;(r) decays 
as 1/r and is consistent with an electrostatic repulsion 
between molecules charged with 0.3 e, as it is described 
in the ab initio results. 

Our study has permitted us to show that a highly- 
ordered chemisorbed phase of single molecules can be 
formed on a metallic surface as a response to the charge 
redistribution upon molecular adsorption. Repulsion 
among localized charges at the molecule/surface region 
are strong enough to hinder nucleation into molecular 
islands, even when H-bonding between molecules is ex- 
pected. We expect similar behavior to occur in other 
molecular systems, thus being a fingerprint of charge 
transfer processes at organic/metal interfaces. In this 
case, the analysis of the nearest-neighbors pair distribu- 
tion of molecules at different coverages turns out to be a 
very valuable tool for the study of adsorbates' interaction 
potential on metallic surfaces. 
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